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ON THE SOLUTION OF A CLASS OF CAUCHY INTEGRAL 

EQUATIONS 


E. DE MICHELI AND G. A. VIANO 


Abstract. In a previous paper we have presented a new method for solving 
a class of Cauchy integral equations. In this work we discuss in detail how 
to manage this method numerically, when only a finite and noisy data set 
is available: particular attention is focused on the question of the numerical 
stability. 


1. Introduction 

In a previous paper we have proved the following theorems. 

Theorem 1. Let us consider the following series 

^ CX3 

— {z = x + iy, x,y gR), (1) 

n—0 

and suppose that the set of numbers {fn}^, fn = (p > 0, n = 0,1; 2,...) 

satisfies the following Hausdorff conditions 

n / \ 2+e 

(n = 0,1,2,. ■•;e>0), (2) 

i=0 

where M is a positive constant, and A is the difference operator defined by: A/„ = 
/„+i - U, A7„ = Ei=oi-'^rOfn+k-m. Then: 

(1) series 0 converges uniformly to a function f{z) analytic in the unit disk 
Dq = {z\ |z| < 1}; 

(2) f{z) admits a holomorphic extension to the “cut-plane” {z G C\ (l,+oo)}; 

(3) the jump function F(x) = —i[f+{x) — f-{x)], (f±{x) = lim,,^o fix ± irj)), 

r )>0 

belongs to L^(l,+oo), and, moreover, if p > 1, it is a function of class 

The e in 0 (missing in Ref. [2) is needed to guarantee the continuity of d(A) 
(Carlsonian interpolation of the a„’s) at A = —1/2 + iu, u G K (see next formula 

(nnji). 

Let us remind the reader that these results in the case p = 0 (i.e. /„ = a„) 
are due to Stein and Wainger j^l; however, these authors do not formulate the 
requirements on the Taylor coefficients a„ by the use of the Hausdorff condition 

0. 

Theorem 2. If in series o the coefficients a-n satisfy the Hausdorff condition 0 
(i.e., fn = (In, n = 0,1, 2,.. .), then the jump function F{x) can be represented by 
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the following expansion, that converges in the sense of -norm: 


OO 


F{x) = ^ Cm4>m{x) {x G [1,+Cxo)), 


( 3 ) 


m—0 

where the coefficients Cm given by: 



( 4 ) 


Pm being the Pollaczek polynomials. The functions 4>m(x) form a basis in L'^{0, +oo), 
and are expressed by: 


0m (^) — Pmijz) 


( 5 ) 


X 


where Bm{x) = i"‘-\/2im(2/x), Lm being the Laguerre polynomials. 

From these results it derives that we can formally solve the following integral 
equation of Cauchy type: 



( 6 ) 


if the infinite set of Taylor coefficients On = f^^\0)/n\ {n = 0,1, 2,...) is known, 
and, in addition, if they are supposed to satisfy the Hausdorff condition with 
fn = ctn- But, at this point, a very serious problem is met: How can this proce¬ 
dure be managed numerically? In fact, it must be noted that, in practice, only 
a finite number of Taylor coefficients can be known, and, moreover, they are usu¬ 
ally affected by numerical errors. Furthermore, if these coefficients are the results 
of experimental measurements the situation is even worse because also the error 
proper of any measurement must be considered. Finally, let us remark that the 
determination of the jump function F{x) starting from an approximate knowledge 
of the function f{z) is a typical example of improperly posed problem in the sense 
of Hadamard: the solution does not depend continuously on the data. In order to 
make this point more clear notice that the problem of solving the Cauchy integral 
equation ® is strictly connected to the problem of the analytic continuation up to 
the cut. We can conformally map the cut z-plane onto the unit disk in the C-plane 
(i.e. the domain |C| < 1); in this map the upper (lower) lip of the cut is mapped in 
the upper (lower) half of the unit circle. Therefore the continuation up to the cut 
corresponds to the continuation up to the unit circle |C| = 1. As is well known, the 
uniqueness of the analytic continuation does not guarantee its continuity in the 
or in the uniform topology. Furthermore, performing the analytic continuation up 
to the boundary of the analyticity domain (i.e. up to the unit circle in the C-plane 
geometry) is a severely ill-posed problem 

In the standard methods of regularization, whose credit is generally due to 
Tikhonov H E], normally some appropriate global bound on the solution are 
imposed, by assuming some prior knowledge on the solution itself. In this way 
a subspace of the solution space is determined, and then a solution belonging to 
this subspace is looked for. If this subspace is compact, then continuity follows 
from compactness. In the numerical analytic continuation in the unit disk (C-plane 
geometry), if one avoids going up to the boundary, this subspace is obtained by 
imposing a global bound at the boundary. This point can be easily understood by 
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observing that from a bound on the functions at |C| = 1 a bound on their deriva¬ 
tives can be derived (in any point inside the unit circle), via the Cauchy integral 
representation. Conversely, if the analytic continuation is required up to the cut 
(i.e. up to the boundary), then a global bound only on the functions is not suf¬ 
ficient, and a uniform bound on the first derivative becomes necessary. Then, in 
this case, the regularization by the use of a-priori bounds is rather cumbersome. 
Furthermore, it must be observed that in several problems of physical interest the 
prior knowledge which allows for imposing bounds on the solution is missing, or can 
be very poor. See on this point Refs. □ini, where the difference between synthesis 
and/or inverse problems is widely discussed: while in the first class of problems the 
prior knowledge is intrinsic to the problem itself, being a part of its formulation, 
for the second class (i.e., inverse problems) this is not true. We are thus motivated 
to go beyond the Tikhonov and Tikhonov-like methods, and, accordingly, to look 
for regularization procedures that do not make use of prior knowledge. 

The method that we present in this paper is very far from the Tikhonov-like 
regularization procedures, and it is actually based on the following observation: 
when in expansion 0 we use the actual data set, which is usually composed by 
only a finite number of data affected by noise, then, although the formal series © 
diverges, nevertheless the effect of the errors (and of the fact that the number of 
data is finite) remains quite small in the beginning of the expansion, and there will 
exist a point (i.e. a certain value of m) where the divergence sets in. Thus, the 
basic idea is to stop the expansion at the point where it turns to diverge. This 
rough and qualitative description can be put in rigorous form by proving that even 
though series © diverges, nevertheless it converges (in the sense of the L^-norm) 
as the number of data tends to infinity and the noise vanishes. This is the main 
point of the paper and it will be proved in Section [2| 

A very delicate point of the method consists in determining a procedure for the 
determination of the truncation point in the approximation which will be proved 
to converge to the unknown jump function F{x). This procedure, that will be 
illustrated in Section|21 which is devoted to the numerical analysis, is based on the 
asymptotic behavior of the Pollaczek polynomial Pm[—i{n + 1/2)] for large values 
of m (at fixed n), whose proof will be given in the Appendix. 


2. Solution of the Cauchy integral equation when the data are a 
FINITE NUMBER OF NOISY TAYLOR COEFFICIENTS 

Let US assume that only a finite number of noisy Taylor coefficients are 
known, and that jon^ — a„| < e (n = 0,1, 2,..., A^; e > 0). Let Cm’^^ be the 
following finite sums 

N 




n—0 


nl 




(7) 


Accordingly, we have Cm°°^ = Cm (see formula ©)■ We can then prove the following 
lemma. 

Lemma 1. The following statements hold true: 

(i) 


El 

m—0 


J0,oo) 




( 8 ) 
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(ii) 


(iii) 


E 

m—O 


2 

= +oo; 


lim = 


N- 


(m = 0,1,2,...); 


(9) 

( 10 ) 


(iv) if kQ{e, N) is defined as 



fco(e, N) = max{fc S N : ^ < C}, 

m—0 

(11) 

then 



lim A:o(e, A^) =+ 00 ; 

(12) 




(v) the sum 


(13) 


m—0 


satisfies the following properties: 

(a) it increases for increasing values of k; 

(b) the following relationships hold true 


M, 




k — 


> 




k —>+00 


(TV!)- 


:{2k) 


2N 


{N fixed). 


(14) 


Proof, (i) In Ref. we proved that the coefhcients a„ in expansion o, which 
are supposed to satisfy condition admit a unique interpolating function a(A) 
(A = tr + iv) which belongs to the Hardy space iJ^(C_i/ 2 ) with norm ||a ||2 = 
sup^>_i/ 2 (E,^ \a(a + iv)\^dv)^^^ (C_i /2 = {A e C | ReA > —1/2}). Furthermore, 
we proved that d{—l/2 + iv), i.e., the restriction of d{a-\-iv) to the line Re A = a = 
— 1/2, belongs to L^(— 00 , + 00 ), and can be represented by the following expansion 
which converges in the sense of the L^-norm: 


a 



00 

^ ^ dm'f^mi^^) ■) 

m—0 


(15) 


where dm = V^Cm (see formula 0), and the functions fimiy) are the Pollaczek 
functions (see [ 3 ]) 

Q + iv) Pm{v), (16) 

(see Refs. Eiiini for what concerns the Pollaczek polynomials). In view of the 
Parseval equality we have: 



(17) 


Next, we recall that a(—1/2 + iv) is the Fourier transform of F{v) exp{v/2), where 
F{v) represents the jump function associated with the cut located at 0 = iv {v € 
(0, + 00 )) in the complex 0-plane geometry (see Fig. lA in Ref. 0 which refers to 
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the series (I/Stt) ™^)- Therefore, in view of the Plancherel theorem and 

equality jni) we get (see formula (46) of Ref. 0): 


^+oo 

d (P + iA 

2 /*+oo 

e^/^F{v) 

/ 

II 

to 

/ —oo 

V 2 J 

^0 



dv 


m—O 


(18) 


Next we pass from the 0-plane geometry of Theorem 2 of Ref. [3| to the z-plane 
geometry of Theorem 2’ of Ref. |2] through the transformation z = exp(—i0). Then, 
to the cut at 9 = iv there corresponds the cut in the z-plane geometry located on 
the real a:-axis {x = Rez) from a: = 1 up to a: = -l-oo. Accordingly, the jump 
function associated with the cut will be F(lna;). Adopting the convention used in 
0 we still denote this latter jump function by F{x) in order to avoid a useless 
proliferation of symbols. Therefore, by the substitution: e" = x, from equality ra 
we obtain: 

-| oo p-\-oc 

\F{^)\^dx = C. (19) 

m=0 

Since from (im equality @ follows. 

(ii) Let us rewrite the sums Cm’^^ as follows: 


AT 




V P 


n=0 




( 20 ) 


where h^n = 1)” ln\. Now, we can write the following inequality: 

1 ' 




N 


V p 


> 


n—0 


^ P 


—i{n + 


-i{ N + 


1 - 




(iV+i)] 


( 21 ) 


In the lemma proved in the Appendix we show that the asymptotic behavior of the 
Pollaczek polynomials Pm[—i{n + 1/2)] for large values of m (at fixed n) is given 
by: 

IM ^ (-1)™^" 


—i n -I- 


-(2m)”. 


( 22 ) 


Therefore, we have: 

Eli.' i-i ("+i)] 


(« + !)] 

Et-„' + 


< 


(iV+i)] 


N-1 

E 

n—0 




(23) 


Wat 


N' 

— (2m) 

Til m—*oo 


\n—N 


^ 0 . 


From (ED, (ED and ED it follows that for m sufficiently large: 




IV! 


(2m)^. 


(24) 
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Therefore, linim^oo = +oo, and statement (ii) follows, 

(ill) We can write the difference as 


N 




^ n—0 


-t \ n-\- - 


Y p 

/ ^ m 


n—N+1 


n\ 




( 25 ) 


In view of the fact that ^ cLnPmi—iin + |)] converges to it 

follows that the second term in bracket (ESJ tends to zero as N ^ oo. Concerning 
the first term we may write the inequality: 


N 




n=0 




N 




n—0 


—i\n + 


(26) 


where the inequalities |a„ — | < e (n = 0,1, 2,..., A^) have been used. Next, by 

rewriting the Pollaczek polynomials + 1/2)] as 


Frr 


—i I n + 


= EpP n + 


J=0 


and substituting this expression in inequality we obtain 


n—0 


m 

(m) 

Ep) 

("4)1 

J=0 



(27) 


(28) 


Next, we perform the limit for N ^ oo. In view of the fact that YJLo pf^\n+l/2y 


is finite, and the series converges, we can exchange the order of the 

sums and write: 


j=0 


n—0 


n - 


(29) 


Finally, performing the limit for e ^ 0, and recalling equality 1(23) , statement (iii) 
is obtained. 

(iv) From dehnition (nj it follows that for ki = fco + 1, we have: 


All 

E| 

m—O 




> c. 


(30) 


Statement (iv) (formula H12ll 'l is proved if we can show that lim^v^cx) fci (e, N) = +oo. 

e —^0 

Let us suppose that lim^v—»oo fci(e, A^) is finite. Then there should exist a hnite 
number K (independent of e and A^), such that limjv —>00 fci(e, N) < K. Then from 
inequality 1301) we have: 


kiie,N) 


m—0 


c< Y. si; 


K 




m—0 


(31) 
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But as N ^ C! 0 ,e —s- 0 we have (recalling also statement (ill) formula dl) 


K 


^<E 


r.(0,Oo) 


m—0 


< 


El 

m—O 


r.(0.Oo) 


= C', 


(32) 


which leads to a contradiction. Then statement (iv) follows. 

(v) Concerning statement (a), it follows obviously from definition 11311 
Finally, the first one of relationships is obvious; the second one follows from 
the asymptotic behavior of Pml—iin + 1/2)] at large m (for fixed n), proved in the 
lemma of the Appendix. □ 


Remark. From statement (v) and formula 113 it follows that, for large values of 
N and small values of e, the sum presents a plateau in a neighborhood of 

k = ko, if the function F{x) being reconstructed is sufficiently regular (see Section 
Eland figs. n3,|33,EK,EK)- 

Now we may introduce the following approximation of the jump function F{x) 
(see expansion ( 0 ): 

ko{e,N) 

F(cW)(a:)= (a; G (l,+oo)). (33) 

m—0 

Approximation F^’^’^^x) is defined through the truncation number ko{e,N); the 
latter can be numerically determined by plotting the sum versus fc, and 

exploiting properties (a) and (b), proved in statement (v) of the previous lemma, 
and the property stated in the remark above. For a more detailed analysis of this 
point see Section El where some numerical examples are given. 

Now, we want to prove that the approximation F^'^'^\x) converges asymptoti¬ 
cally to F{x) in the sense of the L^-norm, as A^ —> oo and e ^ 0. We can prove the 
following theorem. 


Theorem 3. The equality 


lim 

N—^OC 

e-^0 




L2[l,+oo) 


= 0 


holds true. 

Proof. From the Parseval equality it follows that: 




^ m—ko-\-l 


ko 


m—0 


Since 


(0,oo) 

Cm 


= C and limw^oo ko{e, N) = +oo, it follows that 


limw^oo Yl,m=ko + l — 0. 

It is then convenient to rewrite the second term of the r.h.s. of formula (I35II as 
follows. Let us denote by 


€^0 


(34) 


(35) 


( (0,Oo) .r 

(0,00) _ \^rn it m IS even; 

I if TO is odd; 


( 36 ) 
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(gJV) _ 


Cm 


if m is even; 


if m is odd. 


(37) 


Notice that X)m=o - Cm’°°^ = Z)m=o ) and gm’°°^ and 


gm^'^ are real. Next, we introduce the following functions: 

OO 

G(0,oo)(^) ^ ^ g^°--)l[m.m+i[(cr), 


m—0 




(38a) 


(38b) 


m—0 


where I_e is the characteristic function of the set E. From statements (i), (ii), and 
(iii) of the previous lemma (formulae 0, ® and cnj) we obtain 
/• + 00 


= C, 




<ia: = H (aS’”’) 

m—0 
cx> 


G('=’^)(a:) 


N—>oo 
e^O 


m—Q 

(xe [0,+oo)). 


(39) 

(40) 

(41) 


Hereafter we assume, for the sake of simplicity and without loss of generality, that 
every term gm^'^ is different from zero. Next, let X(e, N) be the unique root of the 
equation (x))^ dx = G. Let us indeed observe that (^G^'^’^\x))^ dx 

is a continuous non-decreasing function which is zero for X = 0, and +oo for 
X —> + 00 . Furthermore, from statement (iv) of the previous lemma (formula 11211 1 


we have limAr^oo X(e, N) = +oo. 

e^O, 

Then we can write 


f.X(e.A) ^ 
3 

p+cc 

Jx(e,N) 


G^^'^\x) - (x) 


('g(0,o 


dx 

.A(e.A) 


dx — 2 


G^°’°°\x) fG('’^)(a;) - G^°’°°\x) 


dx. 
(42) 

Next, we perform the limit for N oo and e ^ 0. Concerning the first term at 
the r.h.s. of formula ig2l we have 

p-\-OC 


lim / 

^^■Z^Jx(e,N) 


(G^°'^Hx)y 


dx = 0. 


For what concerns the second term we introduce the following function: 

'g^^’^^x) - G(°’°°^(x) if 0 < a: < X(e,N), 
^0 ifx>X(€,N). 

Then, we have by the use of the Schwarz inequality 




(43) 


(44) 
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Moreover, from dm we have 

-> 0 X € [0,+oo). (46) 

N^oo 

e^O 


The family of functions is bounded in L^[0,+oo), therefore it has a sub¬ 

sequence which is weakly convergent in L^[0, -boo). The limit of this subsequence 
is zero. In fact, let us observe that < 2C; then we consider the func¬ 

tion H^’^’^^{x)^{x), where $(a:) is an arbitrary element of the class of functions 
C'))°(R+). Then, we have |ij('^’'^^(ai)$(x)| < 2C'|$(x)|, and this inequality does not 
depend on N and e. In view of the Lebesgue dominated convergence theorem we 
can then write (see also the limit (gni) 


lim sup 

N —>'00 
e-^0 


^+oo 


H^^'^\x)<i>{x) dx 


= 0 . 


(47) 


Since the set of functions C'))°(M+) is everywhere dense in L^[0,-boo), given an 
arbitrary fnnction ip G L^[0,-boo) and an arbitrary number rj > 0, there exists a 
function S C')j“(R+) such that — ^k\\L'^[o,+oo) < V- Furthermore, through the 
Schwarz inequality we have 


^+oo 


< 


H^'^'^^x) [tpix) - $fe(a;)] 

(•+CO 2^ 


dx 


+ 00 


IV'(t) - $fc(x)| 


1/2 

^' < 2VC'n. 


(48) 


From equality m and inequalities dlHl) we thus conclude that 


lim sup 

N —>00 
e^O 


n-\-oc 


{x)ip(x) dx 


= 0 , 


(49) 


for any ip G L‘^[0, -boo). 

Next, by using the same type of arguments, we can state that if there is an 
arbitrary subsequence belonging to the family that converges weakly in 

L^[0, -boo), then the weak limit of this subsequence is necessarily zero. Finally, from 
the uniqueness of the (weak) limit point it follows that the whole family 
converges weakly to zero in L^[0, -boo). We can thus write 


lim 

N —>oo 
€^0 


r-\-oo 

/ G^°'°°\x)H^^’^\x)dx = 0, 

Jo 


and from inequality (EH) it follows 

rX{e,N) 


lim 

N —>oo 
€^0 ‘ 


G^^'^\x) - G^^'°^\x) 


dx = 0. 


(50) 


(51) 


Since 


771=0 


U,N) _ (O.c 

^771 ^771 


< [Gi'^’^Ha:) - G(°’“)(x)]^ dx, we have 


ko 

lim > 

N —>oo ^ 

e^O 771=0 


rGN) 


_ J0,oo) 



and, in view of equality EH), the theorem is proved. 


(52) 

□ 
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3. Numerical analysis 

In this section various numerical aspects of the issues discussed in the previous 
section will be illustrated. In particular, the effectiveness and the accuracy of the 
reconstruction of the jump function F{x) across the cut will be tested when only 
a finite number of coefficients a„ of the Taylor series o is known. Moreover, tests 
performed with coefficients a„ corrupted by random noise will be illustrated. The 
results hereafter presented summarize a large number of numerical tests performed 
on a variety of sample functions, including smooth functions, regular oscillating 
functions and discontinuous functions. We shall illustrate only the somehow ex¬ 
treme cases, that is, smooth (see Figs. Q] and 0) and discontinuous functions (see 
Fig. 0 , since the behavior of the oscillating ones is very similar to the case of the 
smooth functions. 

In Fig. 0the basic steps for the reconstruction of jump function in the presence 
of only round-off noise are shown (see the legend for the numerical details). In this 
example the test function is Fi{x) = (x € [1, -boo)). This function is 

characterized to be smooth so that we expect the corresponding reconstruction to 
be quite satisfactory. Fig. 04 illustrates a sample of the coefficients a„, which are 
computed as 

/ -l-oo 

x~'^~^Fi(x) dx. (53) 

Then, the coefficients a„ do satisfy condition 0 of Theorem0 In Fig. the plot 
of versus k (see Eq. m) is shown. According to the properties summarized 

in Lemma 0 (see also the remark after the proof of Lemma 0 , a quite extended 
plateau is clearly present before starts increasing rapidly (see Eq. itTUl l. 

This fact allows us to determine the truncation point ko{0,N) which is needed 
for constructing the regularized approximation (see Eq. A simple 

algorithm for the automatic determination of the index k^ has been implemented; 
it exploits both the properties provided by Lemma 0 that fco lies approximately on 
a plateau, and, moreover, that such a plateau must be located before starts 

growing as a power of 2N (see formula d). 

In the general case with e ^ 0, the knowledge of the asymptotic behavior of 
allows for restricting the range of fco by defining an upper limit ka (fco < ka), that 
represents the value of k where approximately the asymptotic behavior sets in; in 
practice, ka is set as the value of k where and its asymptotic behavior starts 

being close enough. Then, the candidate plateaux are located by selecting the 
extended intervals oi k < ka where the modulus of the first numerical derivative of 
is sufficiently small. Einally, k^ is chosen as the largest value of k belonging 
to the interval which is closest, but inferior, to ka- 

Figure 133 shows a quite simple situation, where only one plateau is present in 
the interval from about fc = 15 to fc = 65, and whose value, as expected from Eq. 

is approximately the squared norm of the function Fi (x) (see the horizontal 
dashed line). Panels C and D show two reconstructions of the function Fi(a:); in 
Fig. 0:1 we have fco = 61, which is well within the plateau of Fig. 133 , and the com¬ 
parison between the actual function Fi{x) and its approximation (dotted 

line) shows the good quality of the reconstruction. Similar results (not displayed), 
obtained for different values of ko (15 < kg < 65), indicate that the actual choice 



CAUCHY INTEGRAL EQUATIONS 


11 


B 



C D 




Figure 1: Reconstruction of a smooth jump function. Fi{x) = (x—l)^/(a;'*+a;®+a;^+a;+l), 
X G [l,+oo); the number of coefficients a„ used for the computations is = 30. (A) 
Noiseless coefficients a„ computed according to Eq. I53II . (B) versus k. The 

dashed horizontal line indicates the squared norm of the true function. The algorithm 
for recovering the truncation index gave A:o(0,30) = 61. (C) Reconstruction of the jump 
function. The solid line represents the actual function ^ 1 ( 3 :). The dots are samples of the 
reconstructed jump function computed by using fco(0, 30) = 61. (D) Comparison 

between F\{x) and computed by using A:o(0, 30) = 75. 


of fco is not critical, provided that it is located in the interval corresponding to the 
plateau. Figure CL shows that at fco = 75, that is just after the correct plateau, 
the reconstruction deteriorates evidently. 

Figure [3 illustrates another example of a smooth sample function that gives rise 
to a good reconstruction; in this case ^ 2 ( 2 ;) = (x — 1)^ exp(—x/2) (x € [l,+oo)). 
In Fig. I 3 I) we compare the approximation of f{z) given by /n^z) = 

(Rez G [—1/2,1/2], Imz = 0) with the approximation of f{z) obtained by the 
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A B 



C D 




Figure 2: Reconstruction of a smooth jump function. F 2 {x) = {x — 1)^ exp(—a:/2), 
X € [l,+oo); N = 30. (A) Noiseless coefficients a„. (B) versus k. The dashed 

horizontal line indicates ||F 2 ||^. The algorithm for recovering the truncation index gave 
fco(0, 30) = 82. (C) Comparison between the actual solution F 2 {x) (solid line) and 
(dashed dotted line) computed by using fco(0,30) = 82. (D) Comparison between the 
truncated Taylor series (solid line) and the approximation computed through 

the integral of Cauchy type (see Eq. 0) by using the function (dots), for Re 2 € 

[—1/2,1/2] and Imz = 0. 


Cauchy integral (see Eq. ®) representing the jump function F 2 {x) with its regu¬ 
larized approximation 

For the sake of completeness, it should be mentioned that the erratic behavior of 
the noise can rarely produce very short plateaux located between the true value of 
ko and before starts following its asymptotic behavior (i.e. for k ~ ka)- In 

this case our procedure could fail to recover the correct value of fcg, and this minor 
drawback has been solved heuristically by simply rejecting the plateaux shorter 
than a given threshold L = 5. 
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C D 



Figure 3: Reconstruction of a discontinuous jump function. ^ 3 ( 3 ;) = 1 if a; G [1,10] and 
F 3 {x) = 0 elsewhere; N = 30. (A) versus k. The dashed horizontal line indicates 

jjfsjl^. The algorithm for recovering the truncation index gave A:o(0, 30) = 23. (B) Com¬ 
parison between the actual solution ^ 3 ( 3 ;) (solid line) and (dashed line) computed 

by using fco(0,30) = 23. (C) Comparison between the actual solution F 3 {x) (solid line) 
and (dashed line) computed by using A:o(0, 30) = 44. (D) Comparison between the 

truncated Taylor series En=o (solid line) and the approximation computed through 
the integral of Cauchy type (see Eq. 0) by using the function with fco(0, 30) = 23 

(dashed line), for Re 2 £ [— 1 / 2 , 1 / 2 ] and Imz = 0 . 


However, the situation can become more difficult, and an example is given in 
Fig. 0 Here the jump function ^ 3 ( 3 ;) is a discontinuous one, being constant in the 
range [1,10] and null elsewhere (see Fig. 03). In this case the lack of smoothness 
of F 3 causes the plateau of to be hardly visible, and, consequently, the 

dehnition of the truncation point ko becomes uncertain (see Fig. OK)- In this 
case, our algorithm gave ko = 23, which belongs to the only detectable plateau, 
and the corresponding reconstruction is shown in Fig. 03. Fig. 0F displays the 
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reconstruction of the jump function obtained with ko = 44, that represents the 
index k such that ~ ll^alP (see formula CD). In spite of the fact that, 

according to definition 11111 . fco = 44 is the true truncation index, the corresponding 
approximation is evidently worse than the one obtained with fco = 23 (see Fig. |3j3). 
This fact can be easily explained by recalling that our reconstruction strategy rests 
on a finite truncation of a series expansion (fco < ka < oo for N < oo), and that 
the convergence to the true function is only asymptotic. Then, we pay entirely the 
penalty of the Gibbs-like phenomenon, so that our procedure can lead to inaccurate 
reconstructions in the case of discontinuous jump functions. 

An example of the dependence of the reconstruction on the noise is sketched 
in Fig. 0] Here the test function is again F 2 {x), and the data illustrated in this 
figure are representative of an extensive numerical analysis performed over several 
different test functions that produced similar results. The coefficients a„ have been 
noised by adding white noise, simulated by computer generated random numbers 
uniformly distributed in the interval [— e,e]. Fig. 014 groups the functions 
computed with different values of e in the range from e = 10“^*^ through e = 10“^. 
It is clear that the plateaux become shorter, and, correspondingly, fco decreases 
as e increases. This fact indicates how, when the noise increases, the number 
of coefficients that carry reliable information for reconstructing the jump 

function decreases. Examples of such a reconstructions are plotted in Fig. E^. 
A crude measure of the dependence of the reconstruction stability on the noise e 
is shown in Fig. Efl, where the truncation index fco(e, A^), chosen in this case by 
hand according to its definition CD, is plotted versus the noise bound e. The plot 
indicates that the stability estimate is only logarithmic, as expected for the problem 
of the analytic continuation up to the boundary (see I3IHI). However, in spite of 
this extremely poor stability, nevertheless it is still possible to obtain acceptable 
reconstructions of the jump function, as shown in Fig. EJl. Of course, this is not 
always the case (see for instance the example illustrated in Fig. EJ, since the quality 
of the reconstruction depends strongly on the regularity of the jump function, but, 
at least for some classes of functions, approximation (03) yields useful results even in 
the presence of quite noisy data. This numerical examples make milder some radical 
statements that can be found in the literature (see I3IH!) discouraging any attempt 
of numerical analytic continuation up to the boundary in physical situations, and 
indicate that the numerical procedure here proposed can be successfully applied 
whenever the function being reconstructed is sufficiently regular. To this purpose 
it is worth recalling statement (3) of Theorem ^ In Fig. 0Jl) the plot of the 
reconstruction error, defined as the mean square error of with respect to 

F{x) against the global signal-to-noise ratio (SNR) is shown. It can be seen that 
the reconstruction error is still acceptable for SNR up to about 50 dB, whereas it 
becomes quite large for smaller SNR. 


Appendix 

In this Appendix we prove the following lemma on the asymptotic behavior of 
the Pollaczek polynomials (see also [Bl)- 
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Figure 4: Reconstruction of a smooth jump function from noisy coefficients an\ The 
sample function is F 2 {x); N — 30. The coefficients an have been noised by adding white 
noise uniformly distributed in the interval [—e, e]. (A) versus k computed for e 

ranging from e = 10“^ through e = 10“^° with step 10“^. The rightmost solid line indi¬ 
cates the noiseless (B) Comparison between the actual solution ^ 2 ( 2 ;) (solid line) 

and the reconstructions and F^^^ (dashed lines). (C) Plot of fco(e, 30), set 

by hand according to its definition im . against log(e ^). Each A:o(e, 30) was obtained 
as the average over 300 realizations of the corresponding (D) Plot of the recon¬ 

struction error versus the signal-to-noise ratio SNR. is dehned as the mean square 
error of the reconstruction F^'^’^^x) with respect to the actual function F 2 {x), computed 
in the interval x G [1,20]. The signal-to-noise ratio is defined as the ratio of the mean 
power of the sequence of noiseless coefficients {a„} to the noise variance. 


Lemma 2. The Pollaczek polynomials Pm[—i(n + 1/2)] satisfy the following as¬ 
ymptotic behavior for large values of m (at fixed n), 


Prr 


-i n -b 




■{2my 


( 54 ) 
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Figure 5: Contour used for evaluating integral 


Proof. Let us introduce the following functions Qm^ {y) and (y), that we call 
the associate Pollaczek functions: 

Pm{x)w{x) 


and 


Qm (y) = ^J 
Qin^iy) = 


x-y 
+ Pm{x)w{x) 


■ dx (Imy < 0 ), 

dx (Imy > 0), 


(55) 

(56) 


-oo x-y 

where w{x) = I-r(l/2 + ia;)r(l/2 — ix). Let us now consider the following function 
p[{t, y) defined by 

,_oo X-y 

The series converges to the generating function of the Pollaczek 

polynomials for |t| < 1 |2 E|. The generating function is given by: G{t,x) = 
(1 — Then, we can rewrite integral 15711 as 

(1 - ft)“-i/2(l + ft)-“-V 2 |r(i + ix)r{^ - fa;)} 




"<*•*'*= 2^ 


■ dx 


x-y 


(Imy < 0; |t| < 1). (58) 

We can then integrate the r.h.s. of (EHl) by the contour integration method. Let us, 
indeed, consider the following integral for Ret = 0, 0 < Imt < 1 and Imy < 0, 


Ht^y) = 


i 




^[( 1 - 

zt)(l + ft)] 

t 

_(1 + it)_ 

{-iz-y) 


Gr|i-. 


where 71 is the path shown in Fig. O 


dz, 

( 59 ) 
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Rewriting the term — it)/{\ + it)Y as exp(z In ), it can be easily checked 
that the contributions along the quarters of circle belonging to the path 71 vanishes 
as Rez ^ — 00 . Therefore, through the Cauchy theorem applied to integral ( 1 ^ . 
we obtain, for Ret = 0,0 < Imt < l,Im?/ < 0: 

k 


H{t,y) = 


1 




1 


it + 1 

^ - 1 


(60) 


in view of the fact that the only singularities which play a role are the simple poles 
of r(z + 1/2) at z = —k — 1/2. Let us now recall that 


2 Fi(l,a,a+ l;z) = ^ 


fe=0 


(a + k) 


(|z| < 1, a e C), 


(61) 


where 2^1 denotes the Gauss hypergeometric function; we can thus rewrite the 
expansion as 


Hit,y) = 


1 


2 F 1 ( l,iy + ^Yy + 


3 it + 1 
2’ it-1 


(62) 


(1 - it){iy + 1 / 2 ) 

which shows that H(t,y) is analytic in the half-plane Imt > 0. 

By the use of the following relationship, which holds true for the Gauss hyper¬ 
geometric function 


2 Fi(a, 6 ,c;z) = ( 1 -z) “ 2 -F 1 ( «, c - c; 


z- 1 


z| < 1 , 


z- 1 


we can rewrite the r.h.s. of formula ( 1 ^ as 

1 


H{t,y) = 


F / 1 1 • _L ^ ^ 

2F1 ( l , l , iy +--, —^ 


< 1 , (63) 


(64) 


2{iy+^) 

The hypergeometric series at the r.h.s. of formula (IMIl converges inside a circle of 
radius 2 and center t = i, belonging to the complex t-plane. We can then calcu¬ 
late the following derivatives: H{t,y)/dt'^)t=Q. By the use of the standard 

formula for the derivative of the hypergeometric function we obtain pQ 


1 

m! 


d^H{t,y) 

dt^ 


t=o 


r(m-I-i)r(z?/-I-i) 

2 m+i Y{m -\-iy -\-^) 


3 1 

X 2-Fi ( W -I- l,m -I- l,m-|-iy -I-- 


(65) 


Next, returning to the integral representation of H{t, y) given by formula itKTIl . we 
exchange the integral with the sum. This is legitimate for Ret = 0, |t| < 1, since 
l(f^)“l = I 6 xp[* 2 :ln = 1, (a: € R); then we write, for Ret = 0, |t| < 1 and 

Imi/ < 0: 


if(t,y)= ^t" 
and therefore 


m—0 


QYKy) = - 

ml 


1 

2 

'm Tjl ( 


/*+oo 


Pm{x)w{x) 

x-y 




m—0 


1 (d^H^{t,y) 


dV‘ 


r(TO -I- i)r(zy -I- i) 


X 2 T 1 I w-I-1, TO-|-1, m-I-zy-I- 


z=o Vim + iy+D 

3 r 


2 2 


( 66 ) 


( 67 ) 
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Let US now observe that the function 2 Fi{a,b, c; z)/T{c) is an entire function of 
the variable c; therefore Qm\y) can be analytically continued in the complex y- 
plane, through formula lUTTll and the only singularities that it presents in this (open) 
complex plane are the simple poles of the gamma function T(iy + 1 / 2 ): i.e. y = 
i{N + 1/2), {N = 0,1, 2,...). Now we can employ the following relationship which 
holds true for the hypergeometric function: 

2 Fi{a,b,c]z) = {1- 2 Fi{c- a,c-b,c;z) (|z| < 1), (68) 


and we arrive at 


Q«(2/) 2 F 1 


r(TO + iy + 5) 


11 3 1 

W + + ly + -■ - 


(69) 

Finally, recalling the asymptotic behavior of 2 ^ 1 ( 0 , 6 , c;z) for large |c|, at hxed 
a, b, z (see ^ formula 10, p. 76]), and the asymptotic behavior of the gamma 
function, we obtain 




(70) 


/q\ 

Let us now focus our attention on the function Qm (y), defined for Imy > 0. If 
we change x into —a;, by observing that 


Pm{-x) = {-irPm{x), (71a) 

w{—x) = w{x), (71b) 


we obtain 




/*+oo 


Pm{x)w{x) 
x + y 


dx = (-l)-+iQW(-j/). 


(72) 


This latter relationship, proved for Imy > 0, can be extended by analytic continu¬ 
ation. Next, by the use of the Cauchy formula we obtain 

Pm{x)w{x) 




■dx = -■Kw{jj)Pm{y), 


(73) 


/72 

where 72 is a path which encircles the singularity located at a; = y in a counter¬ 
clockwise sense. From ([721 and ( 1751 ) we finally obtain 

Trw{y)Pmiy) = QmHy) + (- 2 /)- ( 74 ) 

From the asymptotic behavior (TTUI) and by the use of m we obtain: 


Pm{y) 


r(4 +iy)F(i -iy) 

+ (-i)™rQ-*y) i2myy-^/^Y 

Finally, by putting y = —i(n -I- 1/2), we get: 


r ( i + *y ) ( 2 m) 


-* I n -h - 


1 


.r(-n) 

z™(-l) 


(2m)-’^ + (-1)" 


1 


r{n+l) 


{2my 


( 2 m)", 


(75) 


(76) 


which is exactly the result we want to prove. 


□ 
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